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In this paper, we use the cell dynamics method to study the dynamics of phase transformation 
when three phases exist. The system we study is a two-dimensional system. The system is able 
to achieve three phases coexistence, which for simplicity we call crystal, liquid and vapor phases. 
We focus our study on the case when the vapor and crystal phases are stable and can coexist while 
the other intermediate liquid phase is metastable. In this study we examine the most fundamental 
process of the growth of a composite nucleus which consists of a circular core of one phase surrounded 
by a circular layer of second phase embedded in a third phase. We found that there is one special 
configuration that consists of a core stable phase surrounded by another stable phase in a metastable 
liquid environment which becomes stationary and stable. Then, the nucleus does not grow and the 
metastable liquid survives. The macroscopic liquid phase does not disappear even though it is 
thermodynamically metastable. This result seems compatible to the argument of kinetics of phase 
transition developed by Cahn [J. Am. Ceram. Soc. 52, 118 (1969)] based on the construction of a 
common tangent of the free energy curve. 

PACS numbers: 64.60.Qb, 68.55.Ac, 64.75.+g 



I. INTRODUCTION 

The study of the formation of the long-lived metastable 
phase during phase transformation has a long history 0, 
2]. Recently, renewed interest in the metastable phase 
formation, in particular, in the field of soft-condensed 
matter physics 0, Q has emerged because of the rather 
long relaxation time of these materials, which, typically, 
are in the ranges 1 ms to 1 year Although the for- 
mation of the thermodynamically metastable phase is 
not only academically but industrially important because 
many industrial products are in long-lived metastable 
state, the theoretical study of the kinetics of phase trans- 
formation is hindered because of the lack of appropriate 
theoretical and computational models. Hence, a detailed 
understanding of the kinetics of phase transformation us- 
ing realistic modeling is essential. 

The direct microscopic computer simulation of the nu- 
cleation and the kinetics of phase transformation using 
molecular dynamics or the Monte Carlo method is possi- 
ble [3 but is still a difficult task. Even the most funda- 
mental phenomenon like nucleation is still not easy. In 
order to avoid the demand for huge computational re- 
sources, and to get the qualitative (course-grained) pic- 
ture of the kinetics of phase transformation, the meso- 
scopic approach based on the phenomenological model 
called the Cahn-Hilliard Q, Ginzburg-Landau @ or 
phase-field model 0, [ljj, which requires a solution us- 
ing a non-linear partial differential equation, has been 
traditionally employed. Since this approach requires the 
time integration of highly non-linear partial differential 
equations, it is still not easy to simulate the long-time 
behavior of the kinetics of phase transformation [ll| ex- 
cept for the various forms of special analytical traveling 
wave solutions [H GJ, Q E3 . 



In order to understand the full kinetics of phase trans- 
formation with the transient and long-lived metastable 
phase, an efficient simulation method is absolutely nec- 
essary. In this report, we use a formalism which is based 
on the cell dynamics method to investigate the kinetics 
of the metastable phase during the phase transforma- 
tion when the circular grain (nucleus) of the stable phase 
grows. Our result suggests that the metastable phase 
can be long-lived indeed, and it also indicates that the 
cell dynamics method is efficient and flexible enough to 
study the kinetics of the metastable phase during the 
phase transformation. 



II. CELL DYNAMICS METHOD FOR 
THREE-PHASE SYSTEM 

In order to study the phase transformation, it is cus- 
tomary to study the partial differential equation called 
the time-dependent Ginzburg-Landau (TDGL) equation: 
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where ip is the non-conserved order parameter and T 
is the free energy functional (grand potential), which is 
usually written as the square-gradient form: 
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The local part of the free energy h(tp) determines the bulk 
phase diagram and the value of the order parameters in 
equilibrium phases. Traditionally, the double-well form 
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has been used to model the phase transformation of a 
two phase system. 

Puri and Oono transformed this TDGL equation 
l[TJl for the non-conserved order parameter into the space- 
time discretized cell-dynamics equation following a sim- 
ilar transformation of the kinetic equation for the con- 
served order parameter called the Cahn-Hilliard-Cook 
equation [l7| . Their transformation does not correspond 
to the numerical approximation of the original TDGL 
equation. Rather, they aimed at simulating the kinetics 
of phase transformation of real system within the frame- 
work of discrete cellular automata. 

According to their cell dynamics method, the partial 
differential equation Q is transformed into the finite dif- 
ference equation in space and time: 

V(t+l,n)= J F[V(t,n)] (4) 

where the time t is discrete integer and the space is also 
discrete and is expressed by the site index (integer) n. 
The mapping F is given by 

F[ip(t, n)] = -f(ip(t, n)) + [« ${t, n) » -i/>(t, n)] 

(5) 

where f(tp) — dh(ip) / dip and the definition of << * >> 
for the two-dimensional square grid is given by 

i=nn i=nmi 

with "nn" means the nearest neighbors and "nnn" the 
next-nearest neighbors of the square grid. Improved 
forms of this mapping function F for three-dimensional 
case was also obtaine d Il6l ITU Il8| . 

Oono and Puri [l6l Il7| have further approximated the 
derivative of the local free energy f(ip) called "map func- 
tion" by the tanh form: 

/(^) = ^~V-^tanh^ (7) 

dtp 

with A = 1.3, which corresponds to the free energy [l9j : 

h(ip) = -A In (cosh ifj) + ^ip 2 (8) 

and is the approximation to © if A = 1.5 0,0]- Later 
Chakrabarti and Brown |19| argued that this simplifica- 
tion is justifiable since the detailed form of the free 
energy h(ip) is irrelevant to the long-time kinetics and 
the scaling exponent. 

Subsequently, however, several authors used the map 
function f(ip) directly obtained from the free energy h(tp) 
in cell dynamics equation (JSJ as it is |20l |21| and found 
that the cell dynamics equation is still amenable for a 
realistic map function numerically. Ren and Hamley |2lj 
argued that by using the original form of the free en- 
ergy function f(ip) one can easily include the effect of 
asymmetry of free energy and, hence, the asymmetric 



character of two phases can be considered. It is now well 
recognized that this cell dynamics method can reproduce 
the essential feature of the kinetics of phase transforma- 
tion between two phases even though the method is not 
guaranteed [l8j to be an accurate approximation of the 
original TDGL partial differential equation (Q. 

A further extension of the cell dynamics equation to 
the three phase system is simple. One has to introduce 
the free energy function h{tp) of triple-well form, which 
can achieve a three-phase coexistence. In our report, we 
will use one of the simplest analytical forms proposed by 
Widom 

MV0 = ^ + l) 2 (V-l) 2 (^ 2 + e) (9) 

where the parameter e controls the relative stability of 
three phases. Several shapes of the free energy function 
h(ip) for several values of the parameter e are shown in 
Fig-ffl As can be seen from the figure, two phases around 
ip v = —1 which we call vapor and ip c = 1 which we call 
crystal for simplicity always coexist, while another phase 
around tpi = which we call liquid can be metastablc 
when < e < 0.5. The free energy of this metastablc 
liquid phase is higher than that of the stable crystal or 
vapor phases by the amount 

Ah = = 1) - /(^ = 0) = | (10) 

This liquid phase becomes completely unstable and 
disappears when e > e c where e c = 0.5 is the critical 
point. When e < 0, only the liquid phase is stable and 
both the vapor and the crystal phases are metastable 
instead. Since we are interested in the case when only 
one intermediate phase is metastable, we will consider 
the case when < e < e c . As will be shown in the 
next section, even though the metastable liquid minimum 
is irrelevant for the equilibrium phase behavior, it not 
only controls the phase transition kinetics but appears 
as the long-lived macroscopic metastable phase during 
the phase transformation. 

Similar triple-well potentials were used by several 
workers to study the nucleation |2^| and the metastable 
phase formation 0, Q, 0] within the framework of the 
original TDGL or the phase-field model. The importance 
of this triple-well free energy and the appearance of the 
metastable state were also recently suggested experimen- 
tally in a colloid-polymer mixture 0, 01 

III. NUMERICAL RESULTS AND 
DISCUSSIONS 

A. Front velocity of the growing stable phase 

Before looking at the issue of the kinetics of phase 
transformation of a metastable phase in a three-phase 
system, we will briefly look at the growth of one circu- 
lar nucleus from the stable phase after nucleation in a 
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FIG. 1: The model triple-well free energy |23| which can 
achieve a three-phase coexistence. The left phase with ip v — 
— 1 is called vapor, the central phase with ipi = is called liq- 
uid and the right phase with ip c = 1 is called crystal for con- 
venience. The triple point is at e t = 0, when all three phases 
are at equilibrium and can coexist. The liquid phase becomes 
metastable when et < e < e m where e m = 0.2044 though 
the local vapor-liquid or liquid-crystal equilibrium could be 
realized because the common tangent between the stable va- 
por and the metastable liquid phases or those between sta- 
ble crystal and metastable liquid phases can be drawn (bro- 
ken lines). When the parameter e exceeds e m = 0.2044, no 
common tangent can be drawn, so the two-phase equilibrium 
between vapor-liquid and liquid-crystal cannot be achieved. 
When e is larger than e c = 0.5, the liquid phase becomes 
completely unstable. 



two-phase system. In order to simulate the evolution 
of a nucleus, we have to prepare the system as a two- 
phase system in which one phase is stable and another 
is metastable and has higher free energy than the for- 
mer. The free energy difference between the stable and 
metastable phases is controlled by the super saturation in 
usual liquid condensation from vapor and by the under- 
cooling in usual crystal nucleation. Microscopically, this 
free energy difference is necessary for the nucleus of the 
stable phase to grow by overcoming the curvature effect 
of the surface tension [25| . 

In order to study the growth of a stable phase us- 
ing a cell dynamics system, we will consider the time- 
dependent Ginzburg-Landau (TDGL) equation @J of 
square gradient form J2J). The local part of the free en- 
ergy h(ip), which we use is pjj l2l| : 



h(r/>) = i# 2 (l 
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This free energy is shown in Fig. where the liquid phase 
with tpi = is metastable while the crystal phase with 
ip c = 1 is stable, which mimics the liquid-crystal part 
of the triple- well potential !©• The free energy barrier 
between two phases can be tuned by q, while the free 
energy difference Ah between the stable phase at ip c = 1 
and the metastable phase at ipi = can be controlled by 
e. This free energy difference is given by the same formula 
HI U| l as in the model three-phase system. We choose the 
parameter q — 2 in order to mimic the functional form 
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FIG. 2: The model double- well free energy |2Cj, |27|] (solid 
curve) which achieves the two-phase equilibrium is compared 
with the triple- well free energy in Fig. (broken curve) . The 
parameter e determines the free energy difference Ah, and the 
parameter q determines the free energy barrier. We choose 
q — 2 to mimic the triple-well potential. 



The time-dependent Ginzburg-Landau (TDGL) equa- 
tions and J2I for a circular or a spherical growing nu- 
cleus of stable phase in a metastable environment with 
radial coordinate r is written as Il2l 



D 



dr 2 
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dr 



dip 
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where d is the dimension (d = 2 for a circular and d = 3 
for a spherical nucleus) of the problem. Then, the trav- 
eling wave solution having radial symmetry with moving 
interface at R(t) of the form 

^(r,i)=VP0, (13) 
with X = r — R(t) satisfies the differential equation 



d 2 4> ^ # dh 

D dx 2+v dx-d^ = Q 



with 



dR 
~dt 



D(d-1) 
R ' 



(14) 



(15) 



Equation (|14f> represents a mechanical analogue of the 
equation of motion of classical particle in a potential well 
—h subject to a friction force which is proportional to the 
parameter v. Therefore, a finite size of the free energy 
difference Ah is necessary for (I14|) to compensate for the 
dissipation of energy due to friction and to have a solution 
which corresponds to a traveling wave. In other words, 
the moving (growing or shrinking) interface is possible 
only when there is the free energy difference Ah between 
two phases. Therefore, one phase should be metastable 
and another should be stable. 

Equation l|14H has a particular solution only when the 
parameter v takes a specific value. The corresponding 
interfacial velocity dR/dt is given by 



dR 
~dt 



D{d-l) 
R 



(16) 
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where, the second term on the right hand side repre- 
sents the effect of capillary pressure. For a larger nucleus 
with R — > oo, the interfacial velocity becomes v which 
is also the formula for the one-dimension problem with 
d = 1. For a smaller nucleus, the actual interfacial ve- 
locity dR/dt will be smaller than v due to the capillary 
pressure because R > 0. In particular, when dR/dt < 0, 
the nucleus cannot grow. The nucleus with a radius R 
larger than the critical radius R c grows while the one 
with a smaller radius disappears. The critical radius R c 
is determined from dR jdt — which gives 



D(d-1) 



(17) 



Therefore, in two-phase coexistence with v = 0, any cir- 
cular or spherical nucleus with finite radius R disappears, 
and only a flat interface remains. 

The shrinking metastable void within a stable phase 
is also described by the equations |Q to (|17fl . but now, 
with v < and dR/dt < 0. Then, the capillary pres- 
sure in (|l(jfl always accelerates the interfacial velocity, 
and there is no critical radius R c for the void. 

The above steady-state solution of TDGL with a con- 
stant interfacial velocity v was obtained analytically in 
a one-dimension by Chan |12( when the free energy is 
written using the quartic form lllfl . Using his formula, 
the interfacial velocity v of our TDGL model and (J2J) 
with the free energy l|ll|) is given by: 



(18) 



Chan 12] further suggested that if the interfacial width 
is narrow, the interfacial velocity of a circular or spheri- 
cal growing nucleus is asymptotically given by the same 
formula l|18|) . The larger the free energy difference e, and 
the lower the free energy barrier q, the higher the front 
velocity v from (|18fl . 

The critical radius R c of a circular nucleus in a two- 
dimensional system is also given analytically 0, [26] by 




Rc — — 

v 3e 



(19) 



In the metastable environment, a nucleus of a stable 
phase with a radius R smaller than R c shrinks, while 
the nucleus with a radius larger than R c grows and its 
front velocity approaches l|18fl . The void of a metastable 
phase surrounded by a stable environment always shrinks 
regardless of the size of the critical radius. Again, the 
larger the free energy difference e, and the lower the free 
energy barrier q, the smaller the critical radius R c ■ 

We have imported the above free energy (| 1 1 1) into the 
cell dynamics code written by Mathematica TM [2£| for 
the animation of spinodal decomposition developed by 
Gaylord and Nishidate (2^, and simulated the growth 
of the stable crystal phase in a metastable liquid envi- 
ronment and that of a metastable liquid void in a stable 
crystal environment. 



Figures |3] and 0] show the evolution of the stable cir- 
cular crystal phase in a metastable liquid environment, 
and the contraction of a metastable liquid void in a stable 
crystal environment. The system size is 100x100=10000 
and D = 0.5 

The periodic boundary condition is 
used. The initial nucleus or void is prepared by randomly 
selecting the order parameter ip from 0.9 < ip < 1.1 for 
stable crystal and from —0.1 < ip < 0.1 for metastable 
liquid. The initial random distribution is necessary be- 
cause our cell-dynamics system is deterministic and does 
not include random noise. The figures [3] and 0] show that 
the stable circular crystal grows and the metastable cir- 
cular liquid shrinks steadily without changing the circular 
shape appreciably. 





(a) t=0 



(b)t=100 



r i 




(d) t=300 



(c) t=200 



FIG. 3: The evolution of the crystal nucleus (white) in the 
metastable liquid environment (gray) when q = 2 and e = 0.1. 





(a) t=0 



(b)t=100 




(d) t=300 



(c) t=200 



FIG. 4: The contraction of the metastable liquid void (gray) 
in the stable crystal environment (white) when q = 2 and 

6 = 0.1. 
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In Fig.[5]the effective radius r of the circular nucleus of 
a stable crystal grain and a metastable liquid void, which 
is calculated by assuming the circular area S from 




(20) 



is plotted as the function of the time step. We defined 
the area S of crystal as the number of pixels whose order 
parameter ip is larger then 0.5. Figure [S] clearly indicates 
a nearly linear growth of the radius r of the stable phase 
which means the constant front velocity of the liquid- 
crystal interface. 



TABLE I: The front velocities of the growing new stable 
phase (nucleus) and shrinking unstable phase (void) obtained 
from the cell dynamics simulation of Fig. |K| compared with 
the theoretical prediction from ea. 11181 . The critical radii cal- 
culated from eg. (1191 are also shown for reference. 
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(shrinking) 


0.1 


4.71 0.106 


0.023 


0.045 


0.3 


1.57 0.318 


0.133 


0.146 



of the cell dynamics method to the phase transformation 
in a two-phase system and the simulation of the so-called 
Kolmogorov-Johnson-Mehl-Avram (KJMA) kinetic will 
be presented elsewhere p7| . 



._H—a—H J 



Time Steps 

FIG. 5: The evolution and contraction of the effective radii of 
the crystal nuclei and liquid voids calculated from 121 H plotted 
as the function of time steps. □: crystal nucleus when e = 0.1, 
O: crystal nucleus when e — 0.3, A: liquid void when e = 0.1, 
(}: liquid void when e = 0.3. q is set to q = 2 for all four cases. 
The total area is 100x100=10000 pixels. The straight lines 
are the least-square fittings to the numerical data. 



The velocities v estimated from Fig. are summarized 
in Tab. [I] Table[I]shows that the analytical expression in 
ea. (|18fl gives a rough estimate of the front velocity. The 
velocity of shrinking void is always larger than the grow- 
ing nucleus due to the capillary pressure as expected. 
Understanding that the cell dynamics method does not 
attempt to solve the original TDGL directly and, hence, 
is not guaranteed to reproduce the analytical expression 
(|18|l . the discrepancy between the cell dynamics simula- 
tion and theoretical prediction in (|18f) seems not so seri- 
ous. There is also a problem of the definition of the area 
S of the growing phase, which will also numerically affect 
the front velocity calculated from Il20() . 

From the comparison of the cell dynamics simulation 
and the theoretical prediction from the TDGL for the 
two-phase system, we have confidence that this cell dy- 
namics method should be effective to study the qualita- 
tive features of the evolution of the metastable phase in a 
three-phase system, which we will discuss in the next sub- 
section of this paper. More details about the application 



B. Long-lived metastable phase in a three-phase 
system 

Since we are most interested in the evolution or re- 
gression of the metastable phase during the phase trans- 
formation after nucleation, we will study the kinetics 
of phase transformation when three phases exist using 
the model free energy defined by eq.© and depicted in 
Fig. ID 

The phase diagram of this system which we will study 
is shown in Fig. The equilibrium stable vapor phase 
has an order parameter tp v = — 1 and the stable crystal 
phase has ip c = 1, then the vapor-crystal coexistence lines 
(binodal) are the vertical lines at ip v = — 1 and tp c = 1. 
However, the stable liquid at -0; = exists only at the 
triple point e t = 0. 

Once the stable liquid phase disappears from the equi- 
librium phase diagram when e > 0, it can be considered 
to be hidden or buried as a metastable liquid phase within 
the vapor-crystal binodal. Even this metastable liquid 
cannot exist if e > e c where e c = 0.5. The hidden critical 
point is at e c = 0.5 and ip = 0, where the metastable 
liquid becomes completely unstable. 

There are also hidden vapor-liquid and a liquid-crystal 
local coexistence lines (binodals) , which are given by the 
co-tangency points ip v and i\> l of the common tangent 
between the vapor and liquid free energy and by t/'; and 
ip c between the liquid and crystals energy respectively. 
These hidden binodals disappear at e — e m with e m = 
0.2044 where the solutions of the simultaneous equations, 
for example 

h (V) = ti (Vr) 
h (Vr) - h (Vr) = h (V:) - h ^) < (21) 

ceases to exist, and the hidden binodal lines are termi- 
nated by the spinodal lines as shown in Fig. Then, 
the vapor-liquid and the liquid-crystal coexistence could 
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be established locally even when the liquid phase is 
metastable so long as e < e m . 

For the free energy JBJ), the spinodal lines are defined 
by the condition: 



dip 2 



= 



which gives the outer spinodal lines V'spi: 



2 e V7 - 2e + 3e 2 
V ' SP1=± V5~5 + 5^3 



and inner spinodal lines tp sp 2- 



sp2 



V7 - 2e + 3e 2 
5V3 



(22) 



(23) 



(24) 



The latter merge at tpi = and e c = 0.5. Therefore, the 
spinodal region consists of two regions sandwiched by an 
outer spinodal and an inner spinodal lines when e < e c 
while it consists of one region for e > e c as shown in 
Fig. 



vapor-crystal 



stable 
vapor 
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liquid-crystal 



stable 
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FIG. 6: The phase diagram for the model triple-well free 
energy denned by @ . The model realizes the stable vapor and 
crystal phase at %p v = — 1 (solid vertical line) and ip c = 1 (solid 
vertical line) respectively and the metastable liquid phase at 
ipi = (solid vertical line) if < e < e c with e c = 0.5. The 
triple point is at e t = when three phases (■) are all stable 
and can coexist. The metastable liquid ceases to exist at 
ip — and e c = 0.5 which is the hidden critical point (O)- 
There are not only the vapor-crystal binodal (two vertical 
solid lines) but the hidden vapor-(metastable) liquid binodal 
(single-dot chain curve) and hidden (metastable) liquid-crystal 
binodal (double-dots chain curve), which can be constructed 
from the common tangent as shown in Fig.0 These two pars 
of binodal disappear (□) when they terminate the spinodal 
lines (dotted curves) at e m = 0.2044. The spinodal regions 
are also complex because of the existence of the metastable 
liquid phase. 



We have incorporated the above free energy into the 
cell-dynamics code |2j| written by Mathematica TM • 
We have considered the growth of several special forms 
of circular nucleus which consist of two layers of two dif- 
ferent phases embedded in another phase to see the pos- 
sibility of the appearance of the long-lived metastable 
phase. 



In figure we start from the special structure where 
the stable vapor phase is wrapped by the metastable liq- 
uid layer, which is further embedded in a stable crystal. 
The initial crystal, liquid and vapor phases are prepared 
by randomly selecting the order parameter ip from 0.9 to 
1.1 for crystal, from -0.3 to 0.3 for liquid and from -1.1 
to -0.9 for vapor phases. From the analogy of the expan- 
sion of a stable nucleus and the shrinking of metastable 
inner void in a two-phase system, it is expected that the 
outermost crystal phase expands inward and the inner 
vapor core expands outward by consuming the interme- 
diate metastable liquid layer. Actually, the crystal-liquid 
interface move inward as expected while the liquid-vapor 
interface does not move significantly. The vapor core does 
not move significantly even if the radius is larger than the 
critical radius as the surrounding metastable liquid layer 
has finite thickness, while the crystal-liquid front shrinks 
because no critical radius exists. As has been discussed 
in ea.Hlfi[). the capillary pressure accelerates the shrink- 
ing crystal-liquid interface, while it decelerates the grow- 
ing liquid-vapor interface. Then, the slow liquid-vapor 
interface does not have enough time to move because 
metastable liquid is consumed by the fast crystal-liquid 
interface. Finally, the metastable liquid layer disappears 
completely and the stable vapor core is surrounded by the 
stable crystal phase and the vapor-crystal coexistence is 
established. 



VP 



(a) t=0 



(b)t=100 



(d) t=300 



(c) t=200 



FIG. 7: A Gray-level view of the dynamics of a three-layer 
nucleus when e = 0.1. The black area is the vapor phase, 
the gray area is the metastable liquid phase, and the white 
area is the crystal phase, (a) Initial state is the stable vapor 
nucleus (black) wrapped by the metastable liquid layer (gray) 
which is further embedded in the stable crystal phase (white), 
(b), (c) the outer crystal-liquid interface shrinks inward by 
consuming the metastable liquid layer while the inner liquid- 
vapor interface does not expand significantly, (c) Finally the 
stable vapor phase is surrounded by the stable crystal phase 
and the vapor-crystal coexistence is established. 



Figure [S] shows the time evolution of the effective radii 
of the crystal-liquid and the liquid-vapor interface esti- 
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mated from the area of the nucleus from i|20|) . We defined 
the area S of three phases as the number of pixels which 
belong to them; the pixel of vapor is defined by the or- 
der parameter t/j < —0.5, that of crystal by i/j > 0.5 
and that of liquid by —0.5 > tp > 0.5. Figure 03 clearly 
indicates that the effective radius of the crystal-liquid in- 
terface decreases with constant velocity v while that of 
the liquid- vapor interface remains almost constant. The 
crystal-liquid interfacial velocity estimated by fitting the 
straight line to the simulation result is v — 0.026 which 
is again the same order of magnitude as the interfacial 
velocity in a two-phase system listed in Tab.^ 




Time Steps 

FIG. 8: The time evolution of the effective radius of a circular 
crystal-liquid (O) an d liquid-vapor (□) interface when e = 
0.1. The total area is 100x100=10000. The radius of the 
crystal-liquid interface decreases while that of the liquid- vapor 
interface remains almost constant. 



When the metastablc liquid core is surrounded by the 
stable crystal and vapor phases, the crystal-liquid inter- 
face shrinks and the metastable liquid void disappears 
as shown in Fig. [§] The vapor-crystal interface does 
not move appreciably because the vapor and the crys- 
tal phases are in equilibrium and can coexist. Since our 
cellular dynamics model uses continuous order parame- 
ter tp, the order parameter changes continuously in space. 
Then, there is always a thin layer of liquid with ipi = 
between the crystal core with tp c = 1 and the vapor en- 
vironment with ip v = in Fig. 

Figure El shows that the effective radius of the vapor- 
crystal interface remains constant while that of the 
crystal-liquid interface decreases almost linearly. The 
crystal-liquid interfacial velocity is estimated to be v = 
0.023 which is again comparable to the values shown in 
Tab. [I] for the two-phase system. 

These two examples have clearly indicated that the 
behavior of the stable core and the metastable void in 
a triple-phase system is similar to those in a two-phase 
system. 

When the composition of the circular nucleus changes, 
a very interesting behavior is observed. In Figure ITT1 we 
start from the special structure where the stable circular 
crystal core is wrapped by a stable vapor layer, which 
is further embedded in a metastable liquid environment. 
The initial crystal, liquid and vapor phases are prepared 
by randomly selecting the order parameter tp from 0.9 





(a) t=0 



(b)t=100 





(d) t=300 



(c) t=200 



FIG. 9: The same as figure|7|when e = 0.1 but for the differ- 
ent ordering of the layers, (a) Initial state is the metastable 
liquid void (gray) wrapped by the stable crystal (white) layer 
which is further embedded in a stable vapor (black), (b), 
(c) the metastable liquid void shrinks while the outer stable 
crystal-vapor interface remains the same. (c)Finally the sta- 
ble crystal phase is surrounded by the stable vapor phase and 
the vapor-crystal coexistence is established. There is always 
a thin layer of liquid at the solid-vapor interface because we 
use continuous order parameter ip. 



o o o o o o o 
Vapor-Crystal 



O G O 




FIG. 10: The time evolution of the effective radius of circular 
vapor-crystal (O) and crystal-liquid (□) interface when e = 
0.1. The radius of outer crystal- vapor interface remains the 
same while that of the crystal-liquid interface decreases. 



to 1.1 for crystal, from -0.3 to 0.3 for liquid and from 
-1.1 to -0.9 for vapor phases respectively as before. Intu- 
itively, we expect that the stable inner crystal core may 
not grow because of the crystal-liquid coexistence while 
the stable outer vapor layer expands outward by con- 
suming the metastable liquid environment. However, we 
observe that this threefold structure is rather stable for 
a very long time. The crystal phase as well as the vapor 
phase cannot grow and the metastable liquid phase sur- 
vives and occupies almost the same region for a long time 
as if the vapor-liquid coexistence is locally established at 
the liquid-vapor interface. 
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o 



(b) 1=1000 



(a) t=0 



(b) 1=1000 





o 



(d) t=3000 



(c) t=2000 



(d) t=3000 



(c) t=2000 



FIG. 11: The Gray-level view of the dynamics of a three- 
layer nucleus when e = 0.1. (a) Initial state is the stable 
crystal core (white) wrapped by the stable vapor (black) layer 
which is further embedded in the metastable liquid (gray), 
(b), (c), (d) The initial structure does not change appreciably 
even after very long time steps t=1000, 2000, 3000. Note the 
difference of time scale from Figs. |7] and |5] Again, there is 
always a thin layer of liquid at the solid- vapor interface as in 
Fig. H 



Figure IT21 shows the time evolution of the effective ra- 
dius of the vapor-crystal and liquid-vapor interface esti- 
mated from the area of the nucleus calculated from (J2J3J. 
We confirm the stable vapor-crystal as well as the stable 
liquid-vapor interfaces which do not move appreciably. 
In order to confirm the numerical accuracy of Fig. 1111 
we check the evolution of a dual system where crystal 
and vapor is exchanged. Figure ED clearly indicates that 
the dual system shows exactly the same behavior as in 
Fig. ED 



o o o o o 
Liquid-Vapor 



Vapor-Crystal 

□ □□□□□□□□□□□□□ 



Time Steps 

FIG. 12: The time evolution of the effective radii of a circular 
vapor-crystal (□) and liquid-vapor (O) interface when e = 
0.1. The total area is 100x100=10000. The two radii which 
remain almost constant indicating that the growth of stable 
core is prohibited and the long-lived liquid phase survives even 
though it is thermodynamically metastable. Note the long 
time scale compared with Figs. IH1 and 1101 



FIG. 13: The same as Fig. 1111 but now the stable vapor 
core (black) is wrapped by the stable crystal (white) layer 
which is further embedded in the metastable liquid (gray). 
We observe exactly the same morphology as in Fig. 1111 if we 
exchange black (vapor) and white (crystal). 



These striking results may be interpreted from the 
shape of the free energy in Fig. using the argument 
of Cahn Since the parameter e = 0.1, the system re- 
mains in the region of phase diagram (Fig.[SJ) where a hid- 
den liquid-vapor binodal exists and we can draw a com- 
mon tangent between the vapor and liquid phases. This 
means that it is possible to establish local liquid-vapor 
coexistence by changing the local pressure or the chemi- 
cal potential even though the liquid phase is metastable 
(Fig. ID)- Then, the liquid- vapor interface can not move. 
The faceted structure of the liquid-vapor interface ap- 
pears probably due to the high symmetry of the problem 
since we put the nucleus at the center of the area and the 
nearest and the next-nearest neighbors are used to calcu- 
late the Laplacian. The flat interface is also favorable to 
mitigate the capillary pressure which acts to expand the 
liquid-vapor interface. A similar faceted structure ap- 
pears also in the stable liquid- vapor interface in Fig. 

Furthermore, since this stable vapor layer is so tightly 
attracted by an inner crystal core to maintain vapor- 
crystal coexistence, the vapor-crystal interface also can 
not move. Therefore this special three-layer structure 
which Renth et al. 5] called the " boiled-egg crystal" be- 
comes rather stable. The existence of this stable crys- 
tal wrapped by stable vapor layer in metastable liquid 
is predicted theoretically from the shape of the free en- 
ergy El Q an d suggested experimentally in a colloid- 
polymer mixture 0, l30| . 

We note in passing, that this boiled-egg crystal is in 
sharp contrast to the transient metastable phase pre- 
dicted from the steady state solution of the Ginzburg- 
Landau equation. This transient phase appears when it 
is sandwiched by the two stable phase due to the dif- 
ference of the front speed of two interfaces of the stable 
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and metastable phase 0, 0, while our long-lived 
metastable phase appears when it surrounds the two sta- 
ble phases. 

As the parameter e increases further above 0.2, the 
hidden liquid-vapor coexistence cannot be established 
(Fig. m, because we cannot construct a common tan- 
gent. Then this boiled-egg crystal structure cannot re- 
main stable as shown in Fig. 1141 The vapor phase starts 
to grow by consuming the metastable liquid phase, while 
core crystal phase remains almost the same size and the 
shape. Finally the metastable liquid phase disappears 
and the stable crystal phase is wrapped by the stable 
vapor phases and the vapor-crystal coexistence is estab- 
lished. Here again, the initial crystal, liquid and vapor 
phases are prepared by randomly selecting the order pa- 
rameter if) from 0.9 to 1.1 for crystal, from -0.3 to 0.3 for 
liquid and from -1.1 to -0.9 for vapor phases. 




(d) t=300 (c) t=200 



FIG. 14: The same as Fig. [TTlwhen e = 0.3. (a) Initial state 
is the stable crystal core (white) wrapped by the stable vapor 
layer (black) which is further embedded in the metastable 
liquid, (b), (c) The stable vapor layer grows by consuming 
the metastable liquid environment, while the stable crystal 
core remains the same, (c) Finally the stable crystal phase is 
surrounded by the stable vapor phase and the vapor-crystal 
coexistence is established. 

Figure ^] clearly indicates that the effective radius 
of the liquid-vapor interface increases while that of the 
vapor-crystal interface remains the same as the func- 
tion of time. The liquid-vapor interfacial velocity esti- 
mated by fitting the straight line to the simulation data 
is v = 0.112 which is the same order of magnitude listed 
in Tab. U for the two-phase system. 

Similarly, this long-lived boiled-egg structure will be 
destroyed by the thermal noise, which can be simulated 
by using the cell dynamics equation |16( 

ip(t+l,n)=F[ip(t,n)}+Br)(t,n) (25) 

instead of (0J, where B is the amplitude of the noise, 
ry(t, n) is a uniform random number between -1.0 and 
1.0. 




50 100 150 200 250 300 

Time Steps 



FIG. 15: The time evolution of the effective radii of the cir- 
cular vapor-crystal (□) and liquid-vapor (O) interfaces. The 
total area is 100x100=10000. The radius of liquid- vapor in- 
terface increases linearly while that of the vapor-crystal inter- 
face remains almost constant as the function of time. 



In Fig. we start from the same special structure 
as in Fig. 1111 when e = 0.1 but the thermal noise with 
B = 0.07 is included. The thermal noise certainly 
destroys the long-lived metastable configuration as ex- 
pected, but it still remains rather stable for a long time. 
The thermal noise also destroys the circular or even the 
faceted structure and the liquid-vapor interface becomes 
flat. The thermal noise acts to eliminate the curvature 
of the interface to suppress the capillary pressure. 




(d) t=3000 (c) t=2000 



FIG. 16: The same as Fig. llll but the thermal noise with B = 
0.07 is included. The long-lived metastable state is destroyed 
by the thermal noise, however the structure is still rather 
stable for a long time. 



Naturally, the larger the thermal noise, the shorter 
the lifetime of the metastable configuration as shown in 
Fig. El A similar effect of noise on the time-scale of 
the evolution was observed in the TDGL model of nucle- 
ation Q for a non-conserved order parameter. 

These last four examples indicate that the kinetics of 
phase transformation is definitely affected by the pres- 
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FIG. 17: The same as Fig.^Jfor various noise lever B — 0.00 
(Fig.[T2J, B = 0.07, and B = 0.1 when e = 0.1. The effective 
radius of the liquid- vapor interface starts to increases and that 
of the vapor-crystal interface starts to decrease as we increase 
the noise level B. However, the time scale is still much longer 
than Fig. [T51 



ther extension and modification of the method to explain 
the formation of various intermediate phases will be pos- 
sible. 
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ence of metastable phase and the hidden binodals which 
is in no way related to the equilibrium phase diagrams. 
The long-lived metastable phase could appear macro- 
scopically if it accommodates the special composite nu- 
cleus which consists of a stable crystal core surrounded 
by an equally stable vapor layer. 



IV. CONCLUSION 



In this paper, we have used the cell-dynamics method 
to study the evolution of a single composite nucleus. We 
have studied the three-phase system which has a hidden 
binodal with two stable and one metastable phases. We 
have called two stable phases, crystal and vapor, and 
one metastable phase liquid. We could successfully sim- 
ulate the evolution of stable phases and the regression of 
metastable phase. We have found, however, one special 
configuration of a stable crystal core wrapped by a stable 
vapor layer embedded in the metastable liquid environ- 
ment becomes stable and stationary for long time. This 
means that the long-lived metastable environment phase 
can persist and can appear as a macroscopic phase even 
though it is thermodynamically metastable. According 
to the argument of Cahn @], this result can be inter- 
preted from the hidden liquid-vapor binodal which can 
be constructed from the common tangent between the 
stable vapor and metastable liquid phases. 

In conclusion, we have used a cell dynamics method 
to study the growth of a single nucleus which has tra- 
ditionally been explained using the partial differential 
equation derived from time-dependent-Ginzburg-Landau 
equation Q or the so-called phase field model 0, fiof . 
We have found that the long-lived metastable phase can 
appear during the phase transformation as predicted by 
several researchers 0,0, [3, 5J. This cell-dynamics method 
is not only flexible but numerically stable to handle such 
a complex situation when many phases can coexist. Fur- 
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